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Determination of pulsar parallaxes and proper motions addresses fundamental astrophysical open 
issues. Here, after scrutinizing the ATNF Catalog searching for pulsar distances and proper motions, 
we verify that for an ATNF sample of 212 Galactic run away pulsars (RAPs), which currently run 
across the Galaxy at very high speed and undergo large displacements, some gravitational-wave 
(GW) signals produced by such present accelerations appear to be detectable after calibration 
against the Advanced LIGO (LIGO II). Motivated by this insight, we address the issue of the pulsar 
kick at birth, or short rise fling from a supernova explosion, by adapting the theory for emission 
of GW by ultrarelativistic sources to this case in which Lorentz factor is 7 ~ 1. We show that 
during the short rise fling each run away pulsar (RAP) generates a GW signal with characteristic 
amplitude and frequency that makes it detectable by current GW interferometers. For a realistic 
analysis, an efficiency parameter is introduced to quantify the expenditure of the rise fling kinetic 
energy, which is estimated from the linear momentum conservation law applied to the supernova 
explosion that kicks out the pulsar. The remaining energy is supposed to be used to make the star 
to spin. Thus, a comparison with the spin of ATNF pulsars having velocity in the interval 400 — 500 
km s _1 is performed. The resulting difference suggests that other mechanisms (like differential 
rotation, magnetic breaking or magneto-rotational instability) should dissipate part of that energy 
to produce the observed pulsar spin periods. Meanwhile, the kick phenomenon may also occur in 
globular and open star clusters at the formation or disruption of very short period compact binary 
systems wherein abrupt velocity and acceleration similar to those given to RAPs during the short 
rise fling can be imparted to each orbital partner. To better analyzing these cases, pulsar astrometry 
from micro- to nano-arsec scales might be of much help. In case of a supernova, the RAP GW signal 
could be a benchmark for the GW signal from the core collapse. 

PACS numbers: 04.30.Db, 04.80.Nn, 97.60.Gb 



I. ASTROPHYSICAL MOTIVATION 

The detection of gravitational waves (GW) is one of 
the greatest goals of today's relativistic astrophysics. 
Because of the great success of Einstein's general 
relativity (GR) in explaining the dynamical evolution of 
the binary pulsar PSR 1913+16; continuously observed 
over almost 30 years by Taylor and Hulse [ijj, most 
researchers in the field are confident that the very first 
GW signal to be detected must come from coalescing 
binary neutron star-neutron star or black hole-neutron 
star (NS-NS, BH-NS) systems. Ongoing LIGO science 
runs II l| for the first time have set firm upper limits on 



the rate of binary coalescences in our Galaxy: at the level 
of a few events per year. Meanwhile the model-dependent 
estimate of such a rate at cosmic distance scales varies 
around 10 -5 [13, fl3l. [2l| . The expectation is enormous 
as well for an event like the supernova (SN) explosion. 

In this paper we focus on the radiation emitted when 
a pulsar is kicked out from the supernova remnant 
and is accelerated with very high speed, albeit non 
ultrarelativistic. Gravitational radiation is emitted 
whenever the nascent neutron star or pulsar changes its 
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velocity, a condition that is fulfilled by the large thrust 
given to the nascent pulsar during the SN explosion, 
and also due to the pulsar dragging against the SN 
ejecta. (Recall that the gravitational field generated by a 
motion with constant velocity is nonradiative) . We shall 
show that each individual event of launching a nascent 
pulsar to drift or run across the Galaxy, after a SN 
explosion, generates a gravitational wave signal that can 
be detected by current and planned GW interferometers 
such LIGO I, LIGO II, VIRGO, GEO-600, or TAMA- 
300. Nonetheless, it is stressed from the very beginning 
that the strain of the GW signal that we estimate below 
is a lower bound. The reason is that we are not taking 
into account the contribution to the GW amplitude, h, 
stemming from the effective acceleration that the pulsar 
receives at the kick at birth, which could be very large. 
Indeed, the distance over which the RAP is dragged 
against the supernova ejecta could also be astronomically 
large. 

At first glance, it appears that the detection of the 
GW signal from the short rise must have an event rate 
on the order of that for SN events in our galaxy, i. e. ~ 
1 per 30-300 years. In principle such a rate has to be the 
same because the flinging of the pulsar could be just the 
aftermath of most supernova core collapses. However, 
in galactic globular clusters, where the frequency of 
formation and disruption of binary systems harboring 
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compact stars is relatively high, the kick velocity of 
a pulsar partner at the binary disruption could be of 
the same magnitude as the one for galactic RAPs. 
Indeed, Vlemmings et al. fT|| are pursuing a VLBI 
astrometry program to perform measurements of pulsar 
proper motions and parallaxes with spatial resolution of 
microarsecond. As an illustration of the performance 
of the high precision pulsar astrometry system they 
simulated the detection of the disruption of a binary 
system in the stellar cluster Cygnus Superbubble. This 
is a quite estimulating perspective for pulsar astrometry 
in the near future. 

Therefore, there is a potential to detect such a 
disruption or formation of a binary system through the 
GW signal emitted by the abrupt acceleration or sudden 
increase in velocity of a given pulsar. Observations in 
the range of micro- or nano-astrometry with planned 
satellites like GAIA and other instruments on Earth and 
space may help to identify such a phenomenon. An 
advanced technology useful for this purpose is currently 
operative or under planning (see Ref.[l71| and references 
therein). 

The physical reason for advocating in favor of this 
novel possibility for GW detection is based on a number 
of results on the dynamics and kinematics of galactic 
RAPs and also on GR, as is shown next. 

• 1) In a seminal paper, Spruit and Phinney Q 
demonstrated that regardless of the actual physical 
engine driving the natal (instantaneous) pulsar 
kicks, that mechanism can only act upon the pulsar 
over a very restricted time scale, AT~j~, on the 
order of ~ 0.32 s. This is a fundamental piece of 
the dynamics of observed pulsars that seems to be 
overlooked by workers in the field. 

• 2) As shown in the next Section, pulsar surveys 
have shown that the largest (V i max ), smaller 
(K min ) and average (K avo ) spatial (3-D) velocity 
of cataloged pulsars are 5000 km s _1 , 84 km s _1 , 
an d 450 -500 km s _1 , respectively @, H El (see 
Fig 1(c)). That is, V* changes by more than an 



order of magnitude. 

• 3) Observations of NSs in binary systems have 
allowed to estimate their masses with very high 
accuracy. The current values of this physical 
property gather around 1.4 M Q [8|. 

• 4) According to GR the acceleration of a mass 
(M*) to a relativistic velocity (V*); i. e. energy, 
should generate a GW pulse, the nature of which 
has a very distinctive imprint: the so-called 
gravitational wave memory [y, [l8|, [l9| , which makes 
the amplitude of this GW signal (the GW strain) 
not to fade away by the end of the traction phase 
(the kick), but rather leaves it variable in time as 
much as an energy source remains active in the 
source. 



II. GALACTIC RAPS AND SAMPLE 
SELECTION FROM ATNF CATALOG 

Run away pulsars (RAPs) are observed to drift along 
the Galaxy with velocities ranging from 90 km s _1 to 
4000 km s _1 (Australian Telescope National Facility 
(ATNF) Catalog @). Since most stars in our 
galaxy are observed to drift with average velocities of 
15 — 30 km s — it has been suggested that such large 
velocities should be given to the proto-neutron star 
at its birth. The kick velocity imparted to a just- 
born neutron star during the rise time is a fundamental 
piece for understanding the physics of the core collapse 
supernovae, to have an insight on the sources of 
asymmetries during the gravitational collapse and on 
the emission of gravitational waves during this process. 
There seems to be a consensus that the core collapse 
asymmetry driving the kick velocity of a neutron star at 
birth may have origin in a variety of mechanisms (25l . |26| . 

Although the mechanism responsible for the initial 
impulses has not been properly identified yet, a very 
strong constraint on the timescale over which it could 
act on the nascent neutron star (NS) was put forward 
in Ref. 5]. Whichever this thruster engine might be 
it cannot push the star over a timescale longer than 
AT££ - 0.32 s. Combining this with the NS mass, 
which piles-up around 1.4 M Q, and the above average 
velocity, implies that a huge power should be expended 
by the thruster during the pulsar early acceleration 
phase, and hereby a gravitational-wave (GW) burst with 
"memory" (understood as a steady time variation of the 
metric perturbation) must be released during this time, 
according to general relativity [f|. 



The present velocity distribution (see Figsfll 1(c)) is 



constructed based on the analysis of the proper motion 
of radio pulsars in our Galaxy after estimating their 
distances. Pulsar distances are estimated by using annual 
parallax, Hi absorption, or associations with globular 
clusters or supernova remnants, and from dispersion 
measures and galactic electron density models. 

This astronomy field has become a vivid research 
enterprise since the pioneer works of Lyne and Lorimer 
(LL94) @, [|. However, the current view on the 
distribution of kick velocities is far from the original 
one they found through pulsar surveys fifteen years ago. 
LL94 showed that pulsars appear to run around the 
Galaxy (and its halo) with mean spatial (3-D) velocities 
- 450 - 500 km s^ 1 [2], but freeway ("cannon-ball") 
pulsars with velocities as high as 5000 km s" 1 and 
slow-down millisecond pulsars with speeds of 80 km s _1 
wander also around the Galaxy Q. Meanwhile, the 
situation changed dramatically with the pulsar velocity 
survey performed by Arzoumian, Chernoff and Cordes 
(ACC02) [|, who found that such a speed distribution 
appeared to be bimodal indicating one distribution mode 
around 90 km s _1 and another near 500 km s _1 . 
However, this dynamic research field has shown, through 
new studies on pulsar proper motions, that the situation 
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(a) Pulsar number (vertical axis) vs. distance (kpc) to 
Earth (horizontal axis) histogram 




(b) Pulsar number vs. Distance distribution of 212 pulsars 
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(c) Velocity vs. Distance distribution of 212 ATNF 
pulsars having both parameters measured 
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FIG. 1: Data taken from ATNF Pulsar Catalog 



is far from having been settled. For instance, Hobbs et 
al. [l(| have recently challenged the ACC02 findings, 
and have shown that their studies lead to assert that 
there is no evidence of a bimodal distribution. Instead 
they suggest that mean 3-Dim pulsar birth velocity is 
400±40 km s _1 , and that the distribution appears to be 
Maxwellian with 1-Dim rms a = 265 km s _1 , showing 
also a continuum distribution from low to high (present) 
run away speeds, being the last ones the tail of the 
velocity distribution. Hence, looking back on the surveys 
by LL94, the study by Hobbs et al. [1(3] clearly agrees 
with the mean velocity of 450 ±90 km s _1 found by LL94 



A. Sample selection from ATNF database 

In the late times the number of pulsars has increased 
considerably Although various researchers have 

kept update catalogues, since then; in general, these 
have been neither completed nor very accessible. The 
ATNF database is an effort among an australian scientific 
community and some foreing partner institutions. The 
purpose is to compile the largest amount of information 
about pulsars, which is subsequently made available to 
every researcher all over the world. 

The ATNF catalogue [l(| registers all published 
rotation-powered pulsars, including those detected only 
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clicking on the box to the left of the "parameter label" to 
select a couple of pulsar parameters: DIST and VTRANS 
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(a)GW signals from pulsar acceleration and gamma-ray 
burst. Distribution (with the viewing angle) of the GW 
signal produced during the early acceleration phase of a 
pulsar (dashed-green line) with V? vc = 450 km s" 1 , and 
distance = 10 kpc. Comparison with the GW signal from a 
gamma-ray burst (solid-red line), as a function of the jet 
angle to the line-of-sight, with parameters (E) = 10 51 erg, 
7 = 100, and distance = 1 Mpc. 



(see Fig 1(c)). The first parameter is the best estimate 
of the pulsar distance in kpc. The second parameter is 
the transverse velocity, given in km s , based on DIST 
estimates. By default the results will be sorted according 
to the pulsars' J2000 names in ascending alphabetical 
order. However, sorting is possible on any parameter by 
typing the parameter label in the "sort on field text box" 
and selecting whether the sort should be in ascending 
(default) or descending order. However, to obtain as 
much information as possible some careful procedures 
must be followed. Because of this, we selected the 
listing fashion "short format without errors" . The short 
format is used to provide a condensed summary of the 
pulsars parameters (see Table- TV]) . Finally, we selected 
the option "header", that produces a table with header 
information at the top and selected the output option 
"Table" (see Table-fl}. Besides, it is possible to display 
functions of the pulsar parameters as a graph or as a 
histogram. For a normal (x-y) graph, the values to plot 
are defined as regular expressions into the "x, y-axis text 
box" , and the axes of the graph can be displayed linearly 
or logarithmically. The expressions may also contain 
custom-dehned variables or externally defined variables. 



FIG. 2: Angular distribution of the GW signals from pulsar 
acceleration and gamma-ray burst 



at high energies. It also records Anomalous X-ray Pulsars 
AXPs and Soft Gamma-ray Repeaters (SGRs) for which 
coherent pulsations have been detected. However, it 
excludes accretion-powered pulsars such as Her X-l and 
the recently discovered X-ray millisecond pulsars, for 
instance SAX J1808.4-3658 [l|. 

The ATNF catalogue [HJ can be accessed in a number 
of different ways. The simplest one is from a web 
interface [28| allowing listing of the most commonly 
used pulsar parameters, together with the uncertainties 
and information on the references. Several options 
for tabular output format are provided. Currently, a 
total of 67 "predefined parameters" are available. A 
tool is provided for plotting of parameter distributions, 
either as two-dimensional plots or as histograms. Zoom 
tools and interactive identification of plotted points 
are provided. Custom parameters can be defined by 
combining parameters in expressions using mathematical 
operators and functions and these can be either listed 
or plotted. Finally, the sample of pulsars listed 
or plotted can be limited by logical conditions on 
parameters (see Table-|l]), pulsar name (including wild- 
card names) or distance from a nominated position. 
These operational tools are described in more detail 
below and links are provided within the web interface 
to relevant documentation [10( . 

In this particular case, we select from ATNF, by 



III. PRESENT GALACTIC PULSAR 
DYNAMICS HINTS AT DETECTABILITY OF 
GW SIGNALS FROM THE SHORT RISE FLING 
OF RAPS IN SUPERNOVAE 

Next we discuss the prospective for the short rise fling 
of Galactic RAPs in a supernova to generate detectable 
GW signals. The analysis starts from the observed 
kinematic state of RAPs, and goes back to their kicks 
at birth. To such an approach we assume that before the 
kick the star velocity is V® = 0- 

First, let us notice that one can combine the first 
three pieces discussed above together with the extremely 
restrictive character of the thrust timescale (AT"), 
the magnitude of the largest velocity (V J max ), and the 
instantaneous nature of the kick (V® = 0). In so doing, 
one can then conclude that whichever the driving engine 
might be, it should transfer to the star at the kick at birth 
an energy as large as (AE) ~ 10 51 erg. (Recall that 
the NS masses are piled- up around 1.4 M ). Besides, 
keep in mind that according to Newton's second law: 
F* = M^a*, the star's inertia depends directly on its 
rest mass M*. Thus, putting together these ingredients 
one can conclude that the effective acceleration, a*, at 
the kick at birth of each RAP should be the same.[29| 
In other words, because pulsar velocities vary by more 
than an order of magnitude, whilst the inertia (mass) 
of each neutron star is essentially the same (within an 
uncertainty no greater then a factor of 2), one concludes 
that the acceleration at the kick at birth for each of the 
pulsars has definitely to be the same. In all our analysis 
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below we will use this key feature. 

Further, as the kick is an "instantaneous action", 
then the Newtonian kinematics states that the largest 
velocity [3(| should be attained after the longest 
acceleration phase. That is; over the longest traction 
time i^ lck . Therefore, the RAP kinematics during the 
short rise fling follows the law 



* W = K + a * 't 



kick 



(1) 



where V+ n (t) defines the final (at the end of the kick) 
velocity of a given RAP, which is reached after a total 
impulsion time i klck . Consequently, if one assumes for 
the sake of simplicity that the RAP transit through the 
interstellar space does not modify its starting state of 
drifting, then one concludes that the largest (presently 
observed velocity I/ max ~ 5000 km s _1 ) must have 
been reached after the elapsing of the longest timescale 
permitted by the kick mechanism: AT^^ ~ 0.32 s @. 

Thus, after collecting all the pieces on the RAPs 
physics described above one can infer the characteristic 
acceleration received by a typical RAP. To this goal 
on can take the largest present velocity given by Refs. 
B II OS (see Fig.- |l(c)[ ), and the constraint on the 
maximum traction timescale for the kick at birth 
That acceleration then reads 



a* 



^ykick 



10 4 km s" 



(2) 



A. Estimates of present h c , /gw GW characteristics 
hint at prospective detection during RAPs short rise 
fling in supernovae 

The above analysis suggests that the overall time scale 
for acceleration of any given RAP to its present velocity 
(under the assumption stated above) depends exclusively 
on the total energy with which it was thrusted. The 
actual amount of energy expended to thrust the star 
certainly has something to do with the nature of the 
driving mechanism. Hence, the GW dynamics to be 
discussed below, combined with observations of Galactic 
RAPs, may help to enlighten the path to unravel the 
actual thruster. Next we provide order of magnitude 
estimates for both GW amplitude and frequency based 
on parameters of RAPs currently cataloged. (In the next 
Section a rigurous derivation of these GW characteristics 
for the RAP kick at birth is presented). 



must retain information about what happened during its 
short rise fling in the supernova. Thus, the frequency 
of the GW burst emitted during the early impulsion of 
a pulsar in a supernova explosion can be assumed, in a 
very simplified picture, as being inversely related to the 
total acceleration time 18], that is, 



(•kick 

Jew 



^kick 



(3) 



As an example, let us take the values quoted above 
for V™ 3 *, V+ vo , and V\ mm . In these cases, one obtains 
the GW frequencies /g 1 ^ = 3.1 Hz, /g™ = 34.7 Hz, 
and = 186.0 Hz, respectively. It is illustrative to 

compare these frequencies with those provided in Figs.- 
(0 [3 [U [9]) below. As shown next, the last ones would 
be detectable by present-day laser-beam interferometers 
whenever their associate strains (computed below) are 
large enough (see Fig(3] for the estimates of h and /gw 
for the 212 ATNF RAPs sample). 
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FIG. 3: Color Online. LIGO I, LIGO Advanced, VIRGO, 
GEO-600, TAMA-300 strain sensitivities and the h c and f c 
characteristics of the GW signal produced by each of the 212 
pulsars calibrated with respect to LIGO Advanced (green 
squares), and the same sample without calibration (blue 
triangles). The effect of the calibration procedure is evident. 



1. Frequency estimate 

From Eqs.([T])-((2]) one can see that the traction time 
scale for any given pulsar is strictly directly proportional 
to its present spatial velocity. This is a conclusion of 
basic importance. That is, the RAP current kinematics 



2. Amplitude estimate 

To have an order of magnitude estimate of the 
GW amplitude, h 1 from the present kinematic state 
of a galactic RAP one can use the general relativistic 
quadrupole formula 
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„4 r 



(4) 



where r is the pulsar distance as provided by the ATNF 
Catalog, and Ey. = ^M+V+ is the RAP kinetic energy. 

Therefore, if one combines the theoretical prediction 
provided by Eq.([3]) with this approximate GW 
amplitude, one obtains the plot presented in Fig. ([3]). The 
estimates were performed for a sample of 212 pulsares 
selected from the Australian Telescope National Facility 
(ATNF) catalog (U-HU Notice that Eq.® states that 
whenever the RAP velocity is time-varying the signals, 
whose representative points appear above the strain 
sensitivity limit of Advanced LIGO, would be detectable 
by this observatory. For a RAP with no time-varying 
velocity, the signal would be a direct-current-like (DC) 
signal that no detector might observe. (Regarding the 
above estimates, keep in mind that at least the pulsar 
motion around the Galactic center (the bulge) should be 



centripetally accelerated) . 

In Fig.-[3] blue triangles represent GW signals 
computed following this back-of-the-envelope technique. 
Green squares are the obtained after calibrating the 
approximated signals against the strain sensitivity of 
Advanced LIGO. Clearly only a few RAPs would be 
detectable if their GW signals were described in such 
an approximate fashion. Moreover, this graph gives us 
a clear hint at the proper characteristics of the GW 
signal that would be emitted by a run away pulsar 
during its short rise fling in the supernova that makes 
it to drift across the galaxy. These order of magnitude 
estimates suggest that a RAP that was kicked out from 
the supernova ejecta with higher speed and acceleration 
would produce a clearly detectable GW signal. This 
perspective motivated us to study the problem using a 
more rigurous physical description: the theory for the 
emission of gravitational wave signals by ultrarelativistic 
sources, a GW signal which has memory. This is done in 
the next Section. 



TABLE I: Full pulsar sample taken from ATNF. JName is the pulsar 
name based on J2000 coordinates, RAJ right ascension (J2000) in 
hh:mm:ss.s, DECJ decimation (J2000) in +dd:mm:ss, PMRA proper 
motion in right ascension in mas yr" 1 units, PMDEC proper motion in 
declination in mas yr" 1 units, DIST best estimate of the pulsar distance 
in kpc units, VTRANS transverse velocity - based on DIST in km s _1 



# 


PSR[J] 


RAJ 


DECJ 


PMRA 


PMDEC 


DIST VTRANS 




Name 


(hms) 


(dms) 


(mas/yr) (mas/yr) 


(kpc) 


(km/s) 


1 


J0014+4746 


00:14:17.7 


+47:46:33.4 


19.300 


-19.700 


1.82 


237.964 


2 


J0024-7204C 


00:23:50.3 


-72:04:31.4 


5.200 


-3.400 


4.80 


141.384 


3 


J0024-7204D 


00:24:13.8 


-72:04:43.8 


3.800 


-2.400 


4.80 


102.278 


4 


J0024-7204E 


00:24:11.1 


-72:05:20.1 


6.200 


-2.700 


4.80 


153.889 


5 


J0024-7204F 


00:24:03.8 


-72:04:42.8 


4.600 


-2.900 


4.80 


123.747 


6 


J0024-7204G 


00:24:07.9 


-72:04:39.6 


4.200 


-3.700 


4.80 


127.376 


7 


J0024-7204H 


00:24:06.7 


-72:04:06.7 


5.100 


-3.600 


4.80 


142.060 


8 


J0024- 72041 


00:24:07.9 


-72:04:39.6 


4.900 


-3.900 


4.80 


142.515 


9 


J0024-7204J 


00:23:59.4 


-72:03:58.7 


5.350 


-3.420 


4.80 


144.498 


10 


J0024-7204N 


00:24:09.1 


-72:04:28.8 


6.800 


-1.100 


4.80 


156.757 


11 


J0024-7204O 


00:24:04.6 


-72:04:53.7 


4.600 


-2.900 


4.80 


123.747 


12 


J0024-7204U 


00:24:09.8 


-72:03:59.6 


5.100 


-4.500 


4.80 


154.779 


13 


J0030+0451 


00:30:27.4 


+04:51:39.7 


C 


C 


0.30 


8.246 


14 


J0034-0534 


00:34:21.8 


-05:34:36.6 


C 


C 


0.98 


146.280 


15 


J0040+5716 


00:40:32.3 


+57:16:24.9 


c 


c 


4.48 


161.755 


16 


J0117+5914 


01:17:38.6 


+59:14:38.3 


c 


c 


2.14 


186.526 


17 


J0134-2937 


01:34:18.6 


-29:37:16.9 


c 


c 


1.78 


189.641 


18 


J0139+5814 


01:39:19.7 


+58:14:31.7 


-21.000 


-5.000 


2.89 


295.772 


19 


J0147+5922 


01:47:44.6 


+59:22:03.2 


C 


C 


1.91 


96.870 


20 


J0151-0635 


01:51:22.7 


-06:35:02.8 


15.000 


-30.000 


1.93 


306.902 


21 


J0152-1637 


01:52:10.8 


-16:37:52.9 


3.100 


-27.000 


0.79 


101.789 


22 


J0206-4028 


02:06:01.2 


-40:28:04.3 


-10.000 


75.000 


0.88 


315.673 


23 


J0218+4232 


02:18:06.3 


+42:32:17.4 


C 


C 


5.85 


138.673 


24 


J0255-5304 


02:55:56.2 


-53:04:21.3 


0.000 


70.000 


1.15 


381.648 


25 


J0304+1932 


03:04:33.1 


+ 19:32:51.4 


6.000 


-37.000 


0.95 


168.822 


26 


J0323+3944 


03:23:26.6 


+39:44:52.9 


16.000 


-30.000 


1.49 


240.177 


27 


J0332+5434 


03:32:59.3 


+54:34:43.5 


17.000 


-9.500 


1.06 


98.220 


28 


J0335+4555 


03:35:16.6 


+45:55:53.4 


C 


C 


2.08 


63.143 


29 


J0357+5236 


03:57:44.8 


+52:36:57.7 


C 


C 


4.73 


242.561 


30 


J0358+5413 


03:58:53.7 


+54:13:13.7 


9.200 


8.170 


1.10 


64.102 


31 


J0406+6138 


04:06:30.0 


+61:38:40.9 


C 


C 


3.05 


942.673 


32 


J0415+6954 


04:15:55.6 


+69:54:09.8 


c 


c 


1.57 


94.151 



TABLE I: Full sample taken from ATNF Catalog. 



# 


PSR[J] 




RAJ 




DECJ 


PMRA 


PMDEC 


DIST VTRANS 




Name 


(hms) 




(dms) 


(mas/yr) (mas/yr) (kpc) 


(km/s) 


33 


J0437-4715 


04 


:37 


:15 


.8 


-47:15:08.8 


121.310 


-71.530 


0. 


.16 


105.977 


34 


J0452-1759 


04 


:52 


:34 


.0 


-17:59:23.4 


12.000 


18.000 


3 


14 


322.047 


35 


J0454+5543 


04 


:54. 


:07 


.7 


+55:43:41.5 


52.000 


-17.000 





79 


204.903 


36 


J0502+4654 


05 


:02: 


:04 


.5 


+46:54:06.0 


-8.000 


8.000 


1, 


78 


95.475 


37 


J0525+1115 


05: 


:25 


56 


.4 


+ 11:15:19.0 


30.000 


-4.000 


7. 


.68 


1101.985 


38 


J0528+2200 


05 


:28 


:52 


.3 


+22:00:01 


-20.000 


7.000 


2 


.28 


229.047 


39 


J0534+2200 


05: 


:34 


31 


.9 


+22:00:52.0 


-14.700 


2.000 


2. 


.00 


140.669 


40 


J0538+2817 


05 


:38 


:25 


.0 


+28:17:09.3 


-23.530 


52.590 


1. 


.47 


401.682 


41 


J0543+2329 


05 


:43 


:09 


.6 


+23:29:05 


19.000 


12.000 


3. 


.54 


377.151 


42 


J0601-0527 


06 


:01 


:58. 


.9 


-05:27:50.5 


18.000 


-16.000 


7 


,54 


860.899 


43 


J0610-2100 


06 


:10 


:13 


.5 


-21:00:28.0 


7.000 


11.000 


5 


,64 


348.635 


44 


J0612+3721 


06: 


:12: 


:48 


.6 


+37:21:37.3 


C 


C 


1, 


.49 


49.950 


45 


J0613-0200 


06 


:13: 


:43 


.9 


-02:00:47.1 


2.100 


-10.500 





,48 


24.174 


46 


J0614+2229 


06: 


:14 


:17 


.1 


+22:30:36 


-4.000 


-3.000 


4 


,74 


112.361 


47 


J0621+1002 


06: 


:21 


:22 


.1 


+ 10:02:38.7 


3.500 


-0.300 


1. 


.88 


31.310 


48 


J0629+2415 


06: 


:29 


:05 


.7 


+24:15:43.3 


-7.000 


2.000 


4, 


.67 


161.184 


49 


J0630-2834 


06: 


:30 


49 


A 


-28:34:43.1 


-44.600 


20.000 


2 


.15 


498.228 


50 


J0633+1746 


06 


:33 


:54 


.1 


+ 17:46:12.9 


138.000 


97.000 


0. 


.16 


124.953 


51 


J0653+8051 


06: 


:53 


:15 


.0 


+80:52:00.2 


19.000 


-1.000 


3. 


.37 


303.984 


52 


J0659+1414 


06 


:59 


:48 


.1 


+ 14:14:21.5 


44.070 


-2.400 





29 


60.300 


53 


J0700+6418 


07 


:00 


:37 


.8 


+64:18:11.2 


C 


C 


0. 


.48 


22.757 


54 


J0711-6830 


07 


:11: 


:54 


,2 


-68:30:47.4 


-11.000 


19.000 


1. 


.04 


108.249 


55 


J0737-3039A 


07 


:37 


:51 


.2 


-30:39:40.7 


-3.300 


2.600 





,33 


6.639 


56 


J0737-3039B 


07 


:37 


:51 


.2 


-30:39:40.7 


-3.300 


2.600 





,33 


6.639 


57 


J0738-4042 


07 


:38 


:32 


.3 


-40:42:40.9 


-14.000 


13.000 


11.03 


999.054 


58 


J0742-2822 


07 


:42: 


■49 


.0 


-28:22:43.7 


-29.000 


4.000 


1 


,89 


262.313 


59 


J0751+1807 


07 


:51 


:09 


.1 


+ 18:07:38.6 


C 


C 


0. 


.62 


17.809 


60 


J0754+3231 


07 


:54 


:40 


.6 


+32:31:56.2 


-4.000 


7.000 


3. 


92 


149.834 


61 


J0758-1528 


07 


:58 


:29 


.0 


-15:28:08.7 


1.000 


4.000 


3 


.72 


72.717 


62 


J0814+7429 


08 


:14: 


:59 


.5 


+74:29:05.7 


24.020 


-44.000 


0. 


.43 


102.883 


63 


J0820-1350 


08: 


:20 


:26 


.3 


-13:50:55.4 


9.000 


-47.000 


2 


,45 


555.841 


64 


J0823+0159 


08 


:23 


:09 


.7 


+01:59:12.4 


5.000 


-1.000 


1. 


,44 


34.811 


65 


J0826+2637 


08 


:26 


:51 


.3 


+26:37:23.7 


61.000 


-90.000 


0. 


.36 


184.091 


66 


J0835-4510 


08 


:35 


:20 


.6 


-45:10:34.8 


-49.680 


29.900 





,29 


78.542 


67 


J0837+0610 


08 


:37 


:05 


.6 


+06:10:14.5 


2.000 


51.000 


0. 


.72 


174.222 


68 


J0837-4135 


08 


:37 


:21 


.1 


-41:35:14.3 


-2.300 


-18.000 


4. 


.24 


364.772 


69 


J0846-3533 


08 


46 


:06 


.0 


-35:33:40.7 


93.000 


-15.000 


1 


.43 


638.649 


70 


J0908-1739 


09 


:08 


:38 


.1 


-17:39:37.6 


27.000 


-40.000 





,63 


144.143 


71 


J0922+0638 


09 


:22 


14 


.0 


+06:38:23.3 


18.800 


86.400 


1. 


,20 


505.062 


72 


J0943+1631 


09: 


:43: 


:30 


.1 


+ 16:31:37 


23.000 


9.000 


1. 


.76 


206.084 


73 


J0944-1354 


09: 


:44: 


:28. 


.9 


-13:54:41.6 


-1.000 


-22.000 





,69 


72.042 


74 


J0946+0951 


09 


:46 


:07 


.6 


+09:51:55 


-38.000 


-21.000 





,98 


201.720 


75 


J0953+0755 


09: 


:53 


:09 


,3 


+07:55:35.7 


-2.090 


29.460 


0. 


.26 


36.654 


76 


J1012+5307 


10 


:12 


33 


.4 


+53:07:02.5 


2.400 


-25.200 


0. 


.52 


62.407 


77 


J1024-0719 


10 


:24 


:38 


.6 


-07:19:19.1 


-34.900 


-47.000 





,53 


146.072 


78 


J1041-1942 


10: 


:41 


:36 


.1 


-19:42:13.6 


-1.000 


14.000 


3 


,18 


211.606 


79 


J1045-4509 


10: 


:45 


:50 


.1 


-45:09:54.1 


-7.000 


8.000 


3 


,24 


163.287 


80 


J1115+5030 


11 


:15 


:38 


.4 


+50:30:12.2 


22.000 


-51.000 


0. 


.54 


142.196 


81 


J1116-4122 


11 


:16 


43 


.0 


-41:22:43.9 


-1.000 


7.000 


2 


,68 


89.843 


82 


J1136+1551 


11 


:36 


:03 


.2 


+ 15:51:04.4 


-74.000 


368.100 


0. 


.36 


635.733 


83 


J1239+2453 


12 


:39 


:40 


,4 


+24:53:49.2 


-104.500 


49.400 


0. 


.86 


472.410 


84 


J1300+1240 


13 


:00 


:03 


.0 


+ 12:40:56.7 


46.400 


-82.200 


0. 


.77 


344.234 


85 


J1321+8323 


13 


:21 


:46 


.1 


+83:23:38.9 


-53.000 


13.000 


0. 


.77 


199.214 


86 


J1328-4357 


13: 


:28 


:06 


.4 


-43:57:44.1 


3.000 


54.000 


2 


,29 


587.172 


87 


J1430-6623 


14 


:30 


:40 


.8 


-66:23:05.0 


-31.000 


-21.000 


1 


,80 


319.531 


88 


J1453-6413 


14 


:53 


:32 


.7 


-64:13:15.5 


-16.000 


-21.300 


1 


84 


232.391 


89 


J1455-3330 


14 


:55 


:47 


.9 


-33:30:46.3 


5.000 


24.000 





,74 


86.007 


90 


J1456-6843 


14 


:56 


:00 


.1 


-68:43:39.2 


-39.500 


-12.300 





,45 


89.153 


91 


J1509+5531 


15 


:09 


:25 


.6 


+55:31:32.3 


-73.606 


-62.622 


2 


,41 


1104.013 



TABLE I: Full sample taken from ATNF Catalog. 



# 


PSR[J] 




RAJ 




DECJ 


PMRA 


PMDEC 


DIST VTRANS 




Name 


(hms) 




(dms) 


(mas/yr) (mas/yr) (kpc) 


(km/s) 


92 


J1518+4904 


15 


:18 


:16. 


7 


+49:04:34.2 


C 


C 


0. 


.70 


30.914 


93 


J1537+1155 


15: 


:37 


:09. 


9 


+ 11:55:55.5 


1.340 


-25.050 


0. 


.90 


107.038 


94 


J1543-0620 


15: 


:43: 


:30. 


1 


-06:20:45.2 


-17.000 


-4.000 


1. 


15 


95.217 


95 


J1543+0929 


15: 


:43 


:38. 


.8 


+09:29:16.5 


-7.300 


-4.000 


2 


.46 


97.082 


96 


J1559-4438 


15: 


:59 


41 


.5 


-44:38:46.1 


1.000 


14.000 


1 


.58 


105.137 


97 


J1600-3053 


16: 


:00 


:51. 


.9 


-30:53:49.3 


-1.020 


-6.700 


2 


.67 


85.788 


98 


J1603-7202 


16: 


:03 


35 


6 


-72:02:32.6 


-2.800 


-6.000 


1 


.64 


51.481 


99 


J1604-4909 


16: 


:04 


:22. 


9 


-49:09:58.3 


-30.000 


-1.000 


3 


.59 


510.886 


100 


J1607-0032 


16: 


:07 


:12. 


1 


+00:32:40.8 


-1.000 


-7.000 





.59 


19.779 


101 


J1623-2631 


16 


:23: 


:38. 


2 


-26:31:53.7 


-13.400 


-25.000 


2 


,20 


295.848 


102 


J1640+2224 


16 


:40 


:16. 


7 


+22:24:08.9 


1.660 


-11.300 


1. 


19 


64.436 


103 


J1643-1224 


16: 


:43 


:38. 


.1 


-12:24:58.7 


3.000 


-8.000 


4 


.86 


196.863 


104 


J1645-0317 


16: 


:45 


:02. 





-03:17:58.3 


-3.700 


30.000 


2 


.91 


417.022 


105 


J1709-1640 


17 


:09 


:26. 


.4 


-16:40:57.7 


3.000 


0.000 


1. 


.27 


18.063 


106 


J1709+2313 


17 


:09 


:05. 


.7 


+23:13:27.8 


-3.200 


-9.700 


1 


.83 


88.618 


107 


J1713+0747 


17 


:13 


49 


.5 


+07:47:37.5 


4.917 


-3.933 


1. 


.12 


33.541 


108 


J1720-0212 


17 


:20 


:57. 


2 


-02:12:23.9 


-1.000 


-26.000 


5. 


.41 


667.357 


109 


J1722-3207 


17 


:22 


:02. 


9 


-32:07:45.3 


-1.000 


-40.000 


3 


.18 


603.239 


110 


J1735-0724 


17 


:35 


:04 


9 


-07:24:52.4 


-2.400 


28.000 


4. 


.32 


575.570 


111 


J1738+0333 


17 


:38 


:53. 


.9 


+03:33:10.8 


5.600 


4.000 


1. 


.97 


64.275 


112 


J1740+1311 


17 


:40 


:07. 


3 


+ 13:11:56.6 


-22.000 


-20.000 


4 


,77 


672.375 


113 


J1741-3927 


17 


:41 


:18. 


.0 


-39:27:38.0 


20.000 


-6.000 


4 


.75 


470.223 


114 


J1744-1134 


17 


:44 


:29. 


.4 


-11:34:54.6 


19.600 


-7.000 





,48 


46.986 


115 


J1745-3040 


17 


:45 


56 


3 


-30:40:23.5 


6.000 


4.000 


2 


,08 


71.110 


116 


J1752-2806 


17 


:52: 


:58. 


6 


-28:06:37.3 


-4.000 


-5.000 


1 


,53 


46.446 


117 


J1801-2451 


18: 


:00 


59 


.8 


-24:51:53 


2.000 


-3.000 


4 


,61 


78.802 


118 


J1801-1417 


18 


:01 


:51. 


.0 


-14:17:34.5 


-8.000 


-23.000 


1 


,80 


207.810 


119 


J1803-2137 


18 


:03 


:51. 


.4 


-21:37:07.3 


11.600 


14.800 


3 


.94 


351.252 


120 


J1807-0847 


18: 


.07 


:38. 





-08:47:43.2 


-5.000 


1.000 


3 


.60 


87.027 


121 


J1820-0427 


18: 


:20 


:52. 


6 


-04:27:38.1 


C 


C 


2 


,45 


164.266 


122 


J1823+0550 


18: 


:23: 


:30. 


.9 


+05:50:24.3 


5.000 


-2.000 


3 


01 


76.848 


123 


J1824-1945 


18 


:24 


:00. 


4 


-19:45:51.7 


-12.000 


-100.000 


5 


,20 


2482.989 


124 


J1824-2452 


18 


:24 


:32. 





-24:52:11.1 


-0.900 


-4.600 


4. 


,90 


108.887 


125 


J1825-0935 


18: 


:25 


:30. 


.5 


-09:35:22.1 


-13.000 


-9.000 


1 


,00 


74.961 


126 


J1829-1751 


18 


:29 


:43. 


1 


-17:51:03.9 


22.000 


-150.000 


5 


,49 


3945.954 


127 


J1832-0827 


18 


:32: 


:37. 


.0 


-08:27:03.6 


-4.000 


20.000 


4. 


,75 


459.311 


128 


J1833-0827 


18: 


:33 


:40. 


.3 


-08:27:31.2 


C 


C 


5 


,67 


901.625 


129 


J1835-1106 


18 


:35 


:18. 


.2 


-11:06:15.1 


27.000 


56.000 


3 


,08 


907.804 


130 


J1836-0436 


18: 


:36 


:51. 


.7 


-04:36:37.6 


C 


C 


4 


,62 


265.780 


131 


J1836-1008 


18: 


:36 


:53. 


.9 


-10:08:08.3 


18.000 


12.000 


5 


39 


552.813 


132 


J1840+5640 


18 


:40 


I | 


6 


+56:40:55.4 


-30.000 


-21.000 


1. 


,70 


295.141 


133 


J1841-0425 


18 


:41 


:05. 


.6 


-04:25:19.6 


C 


C 


5 


17 


285.843 


134 


J1844+1454 


18: 


:44 


54 


.8 


+ 14:54:14.1 


-9.000 


45.000 


2. 


,23 


485.178 


135 


J1850+1335 


18 


:50 


:35. 


.4 


+ 13:35:58.3 


C 


C 


3. 


,14 


240.040 


136 


J1857+0943 


18 


:57 


:36. 


3 


+09:43:17.2 


-2.940 


-5.410 





,91 


26.537 


137 


J1900-2600 


19: 


:00 


:47 


5 


-26:00:43.8 


-19.900 


-47.300 


2 


,00 


486.568 


138 


J1902+0615 


19 


:02: 


:50. 


2 


+06:16:33.4 


c 


C 


9. 


63 


368.086 


139 


J1905-0056 


19: 


:05 


:27. 


.7 


-00:56:40.9 


C 


C 


6. 


,91 


296.655 


140 


J1906+0641 


19: 


:06 


:35. 


2 


+06:41:02.9 


C 


C 


9. 


.23 


269.607 


141 


J1907+4002 


19: 


:07 


:34. 


.6 


+40:02:05.7 


11.000 


11.000 


1. 


.76 


129.804 


142 


J1909+1102 


19 


:09 


:48. 


6 


+ 11:02:03.3 


c 


c 


4. 


31 


192.770 


143 


J1909-3744 


19 


:09 


:47. 


.4 


-37:44:14.3 


-9.490 


-36.060 


1. 


,14 


200.885 


144 


J1910-5959A 


19 


:11 


12 


7 


-59:58:26.9 


-3.300 


-3.600 


4. 


00 


92.613 


145 


J1910-5959C 


19 


:11 


:05. 


.5 


-60:00:59.7 


-4.100 


-4.600 


4 


,00 


116.855 


146 


J1911-1114 


19 


:11 


:49. 


2 


-11:14:22.3 


-6.000 


-23.000 


1 


,59 


179.179 


147 


J1913+1400 


19 


:13 


2 1 


3 


+ 14:00:52.7 


C 


C 


5. 


.07 


326.934 


148 


J1913-0440 


19: 


:13 


:54. 


1 


-04:40:47.6 


7.000 


-5.000 


3 


.22 


131.322 


149 


J1915+1009 


19: 


:15 


:29. 


.9 


+ 10:09:43.7 


C 


C 


5. 


.32 


177.003 


150 


J1915+1606 


19 


:15 


:27. 


.9 


+ 16:06:27.4 


-2.560 


0.490 


7. 


.13 


88.107 



TABLE I: Full sample taken from ATNF Catalog. 



# 


PSR[J] 




RAJ 




DECJ 


PMRA 


PMDEC 


DIST VTRANS 




Name 


(hms) 




(dms) 


(mas/yr) (mas/yr) (kpc) 


(km/s) 


151 


J1916+0951 


19: 


:16 


:32. 


.3 


+09:51:25.9 


C 


C 


2.88 


123.257 


152 


J1917+1353 


19: 


:17 


:39. 


7 


+ 13:53:56.9 


0.000 


-6.000 


4.07 


115.774 


153 


J1919+0021 


19: 


:19 


:50. 


6 


+00:21:39.8 


-2.000 


-1.000 


3.32 


35.196 


154 


J1921+2153 


19: 


:21: 


1 I 


.8 


+21:53:02.2 


17.000 


32.000 


0.66 


113.382 


155 


J1926+1648 


19: 


:26 


45 


3 


+ 16:48:32.7 


C 


C 


7.74 


597.353 


156 


J1932+1059 


19: 


:32: 


13. 


.9 


+ 10:59:32.4 


94.090 


42.990 


0.36 


177.051 


157 


J1935+1616 


19: 


:35 


:47. 


.8 


+ 16:16:40.2 


-1.000 


-13.000 


7.93 


490.190 


158 


J1937+2544 


19: 


:37 


:01. 


2 


+25:44:13.6 


c 


C 


2.76 


232.973 


159 


J1939+2134 


19: 


:39 


:38. 


.5 


+21:34:59.1 


-0.460 


-0.660 


3.60 


13.731 


160 


J1941-2602 


19: 


:41 


:00. 


.4 


-26:02:05.7 


12.000 


-10.000 


4.74 


351.027 


161 


J1944+0907 


19: 


:44: 


:09. 


3 


+09:07:23.2 


12.000 


-18.000 


1.28 


131.280 


162 


J1946-2913 


19: 


:46 


:51. 


7 


-29:13:47.1 


19.000 


-33.000 


4.31 


778.087 


163 


J1946+1805 


19: 


:46 


53 





+ 18:05:41.2 


1.000 


-9.000 


0.85 


36.492 


164 


J1948+3540 


19: 


:48 


:25. 


.0 


+35:40:11.0 


-12.600 


0.700 


7.87 


470.849 


165 


J1952+3252 


19: 


:52 


:58. 


.2 


+32:52:40.5 


-24.000 


-8.000 


2.50 


299.845 


166 


J1954+2923 


19: 


:54: 


:22. 


.5 


+29:23:17.2 


25.000 


-36.000 


0.42 


87.273 


167 


J1955+2908 


19: 


:55 


:27. 


.8 


+29:08:43.5 


-1.000 


-3.700 


5.39 


97.941 


168 


J1955+5059 


19: 


:55 


:18. 


7 


+50:59:55.2 


-23.000 


54.000 


1.80 


500.880 


169 


J1959+2048 


19: 


:59 


:36. 


.7 


+20:48:15.1 


-16.000 


-25.800 


1.53 


220.211 


170 


J2002+4050 


20 


:02: 


I | 





+40:50:53.9 


C 


C 


8.30 


524.995 


171 


J2004+3137 


20 


:04: 


:52. 


2 


+31:37:10.0 


C 


C 


8.94 


473.870 


172 


J2013+3845 


20 


:13 


:10. 


3 


+38:45:43.3 


-32.100 


-25.000 


13.07 


2521.130 


173 


J2018+2839 


20 


:18 


:03. 


.8 


+28:39:54.2 


-2.600 


-6.200 


0.97 


30.945 


174 


J2019+2425 


20 


:19 


31. 


9 


+24:25:15.3 


-9.410 


-20.600 


0.91 


97.708 


175 


J2022+2854 


20 


:22 


:37. 





+28:54:23.1 


-4.400 


-23.600 


2.70 


307.605 


176 


J2022+5154 


20 


:22: 


:49. 


.8 


+51:54:50.2 


-5.230 


11.500 


2.00 


119.788 


177 


J2023+5037 


20 


:23 


II 


.9 


+50:37:34.8 


C 


C 


1.80 


187.936 


178 


J2046-0421 


20 


:46 


:00. 


.1 


-04:21:26.0 


9.000 


-7.000 


3.83 


207.032 


179 


J2046+1540 


20 


:46 


:39. 


3 


+ 15:40:33.6 


-13.000 


3.000 


2.56 


161.926 


180 


J2048-1616 


20 


:48 


35 


.4 


-16:16:43.0 


117.000 


-5.000 


0.64 


355.328 


181 


J2051-0827 


20 


:51 


:07. 


.5 


-08:27:37.7 


5.300 


0.300 


1.28 


32.214 


182 


J2055+2209 


20 


:55 


:39. 


.1 


+22:09:27.2 


c 


C 


2.15 


52.413 


183 


J2055+3630 


20 


:55: 


31. 


3 


+36:30:21.4 


-3.000 


3.000 


5.56 


111.835 


184 


J2108+4441 


21: 


:08 


:20. 


4 


+44:41:48.8 


3.500 


1.400 


5.28 


94.362 


185 


J2113+2754 


21: 


:13 


:04. 


3 


+27:54:02.2 


-23.000 


-54.000 


1.39 


386.791 


186 


J2113+4644 


21: 


:13 


:24. 


3 


+46:44:08.7 


c 


C 


4.99 


308.455 


187 


J2116+1414 


21: 


:16 


13. 


7 


+ 14:14:21.0 


8.000 


-11.000 


4.43 


285.665 


188 


J2124-3358 


21: 


:24 


:43. 


.8 


-33:58:44.6 


-14.400 


-50.000 


0.25 


61.670 


189 


J2129+1210A 


21: 


:29 


:58. 


2 


+ 12:10:01.2 


-0.260 


-4.400 


8.90 


185.980 


190 


J2129+1210B 


21: 


:29 


:58. 


6 


+ 12:10:00.3 


1.700 


-1.900 


8.90 


107.576 


191 


J2129+1210C 


21: 


:30 


:01. 


2 


+ 12:10:38.2 


-1.300 


-3.300 


10.00 


168.154 


192 


J2129-5721 


21: 


:29 


:22. 


7 


-57:21:14.1 


7.000 


-4.000 


2.55 


97.468 


193 


J2145-0750 


21 


:45 


:50. 


4 


-07:50:18.3 


C 


C 


0.50 


33.461 


194 


J2149+6329 


21 


:49 


:58. 


.5 


+63:29:43.5 


14.000 


10.000 


13.65 


1113.385 


195 


J2150+5247 


21: 


:50 


:37. 


7 


+52:47:49.6 


C 


C 


5.67 


231.241 


196 


J2157+4017 


21: 


:57 


:01. 


.8 


+40:17:45.8 


17.800 


2.800 


5.58 


476.682 


197 


J2219+4754 


22. 


:19 


:48 


.1 


+47:54:53.9 


-12.000 


-30.000 


2.45 


375.304 


198 


J2225+6535 


22 


:25 


:52. 


4 


+65:35:34.0 


144.000 


112.000 


2.00 


1729.770 


199 


J2229+6205 


22 


:29 


II 


.8 


+62:05:36.0 


C 


C 


5.70 


270.235 


200 


J2229+2643 


22 


:29 


:50. 


.8 


+26:43:57.7 


1.000 


-17.000 


1.43 


115.452 


201 


J2235+1506 


22 


:35 


:43. 


7 


+ 15:06:49.0 


15.000 


10.000 


1.15 


98.289 


202 


J2257+5909 


22 


:57 


:57. 


.7 


+59:09:14.8 


C 


C 


6.40 


540.227 


203 


J2257+5909 


22 


:57 


:57. 


.7 


+59:09:14.8 


C 


C 


6.40 


540.227 


204 


J2305+3100 


23: 


:05 


:58. 


.3 


+31:00:01.7 


2.000 


-20.000 


3.92 


373.546 


205 


J2308+5547 


23: 


:08 


13. 


.8 


+55:47:36.0 


-15.000 


0.000 


2.42 


172.097 


206 


J2313+4253 


23: 


:13 


:08. 


.5 


+42:53:12.9 


C 


C 


0.95 


99.296 


207 


J2317+1439 


23 


:17 


:09. 


.2 


+ 14:39:31.2 


-1.700 


7.400 


1.89 


68.034 


208 


J2321+6024 


23: 


:21 


:55. 


2 


+60:24:30.7 


-17.000 


-7.000 


3.21 


279.789 


209 


J2322+2057 


23: 


:22 


:22. 


.3 


+20:57:02.9 


-17.000 


-18.000 


0.78 


91.557 
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TABLE I: Full sample taken from ATNF Catalog. 



# 


PSR[J] 


RAJ DECJ 


PMRA 


PMDEC 


DIST VTRANS 




Name 


(hms) (dms) 


(mas/yr) 


(mas/yr) 


(kpc) 


(km/s) 


210 


J2326+6113 


23:26:58.6 +61:13:36.4 


C 


C 


4.82 


433.576 


211 


J2330-2005 


23:30:26.8 -20:05:29.6 


74.700 


5.000 


0.49 


173.922 


212 


J2337+6151 


23:37:05.7 +61:51:01.6 


C 


C 


2.47 


124.481 


213 


J2354+6155 


23:54:04.7 +61:55:46.7 


22.000 


6.000 


3.31 


357.846 



From the analysis of this Table one can realize that 
most of the SNR that were obtained gather about an 
S/N ~ 0(1). However, an SNR of the order of unity 
is not enough for detection. The statistics which is the 
output of the matched filter must stand over the GW 
detector noise. Thus, if one assumes Gaussian noise in 
the detector and a Gaussian distributed statistics, then 
an SNR of something like 5, what means a chance of 
one in a million, is required to state that it is signal and 
not noise. Notwithstanding, RAPs current kinematics 
suggest that the dynamical situation could have been 
dramatically different during the short rise fling that 
launches each pulsar to drift across the Galaxy. In virtue 
of the potentiality of this phenomenon for the emission 
of powerful gravitational radiation signals, in the next 
Section we analyze the pulsar kick at birth within the 
framework of a theory that appears more appropriate to 
describe that process, and reestimate the characteristics 
of the GW released. 



IV. GRAVITATIONAL WAVES FROM PULSAR 
ACCELERATION 

Observations of extragalactic astrophysical sources 
suggest that blobs of matter are ejected with 
ultrarelativistic speeds (Lorentz factor 7 >> 1) in 
various powerful phenomena involving quasars, active 
galactic nuclei, radio-galaxies, supernova explosions, 
and microquasars. For these sources the emission of 
gravitational waves when such an ultrarelativistic blob is 
ejected from the core of its host source have been studied 
by Segalis and Ori[l9[ (see also for the analysis of 
similar phenomena in gamma-ray bursts). 

Supernova explosions are known to be the natural 
place of birth of neutron stars and pulsars, e. g. The 
Crab Nebula. In processes such these, the nascent pulsar 
is kicked out from the supernova remnant with speed that 
is not so large as compared to that for the blobs in AGN 
or quasars, where a significant fraction of the speed of 
light is typical for the ejecta. Nonetheless, the velocity 
at birth can be considered large if one bears in mind 
that the typical mass of neutron stars is about one solar 
mass [8]. Although a galactic RAP cannot be considered 
an ultrarelativistic source, the pulsar rise fling during the 
supernova can kick it out with very high speed. That is 
why in the discussion below we analyze the gravitational 
radiation emitted during the kick at birth of RAPs within 
the framework of Ref . [19] . Besides, it is also of worth 



to keep in mind that such an approach can be applied 
to the early RAPs dynamics because in this case the 
limit 7 ~ 1 applies. This limit corresponds to the GW 
strain h provided by the well-known general relativistic 
quadrupole formula for the same dynamical situation [|| . 
But bear in mind that this does not describe properly 
the GW angular distribution. 

The pulsar is envisioned hereafter as a "particle" of 
mass moving along the worldline r (t) (with t the 
proper time) and having an energy-stress tensor 




where V a — dr a jdr is the particle 4- vector 
velocity. (Upper[sub])scripts are raised [lowered] with the 
Minkowski metric n^ v ). 

Since we are using the linearized Einstein's equations 
to describe the emission of GW during the early 
acceleration phase of any pulsar, then the resulting 
metric perturbation produced by the RAP can 
straightforwardly be evaluated at the retarded time, 
which corresponds to the intersection time of r a (r) with 
the observer's past light-cone, and then be transformed 
to the Lorentz gauge where it reads: h^ v — h^ — ^n^Jh^. 
This strain can equivalently be written as 

h ^= -vX-rW + ■ (6) 

Notice that is this factor — V\ ■ [x — r(r)] A ; which 
depends on the velocity in the denominator of Eq.©, 
that is responsible for the nonvanishing amplitude (GW 
"memory" 6]) of the GW signal produced by the 
launching of the pulsar into its present trajectory. 

Eq.® must be finally rewritten in the tranverse- 
traceless (TT) gauge, i. e. — > /iTJ; which is the 
best suited to discuss the GW detector's response to that 
signal. This procedure leads to the result presented in 
Eq.Q. 

A detailed analysis (see Refs.p, LL9|) shows that the 
maximum GW strain in the detector is obtained for a 
wavevector, n, orthogonal to the detector's arm, in the 
v - ft plane. Here ft is the unit spatial direction vector 
from the (retarded) pulsar position to the observer, and 
v the unit spatial direction vector defining the 3-D pulsar 
velocity V. In this case, the GW amplitude generated by 
the pulsar kick becomes (c.g.s. units recovered) 
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<(*) 



7(i)Af*/? 2 (t) / sin 3 9 sin 2 A</> 



AQ(1 - (3 cos 9) 



de 



(7) 



allows us to neglect the radial component of the pulsar 
velocity because it does not contribute to the effective 
GW amplitude of the signal being emitted by the RAP, 
or in other words, to the signal to be detected on Earth. 



f c 

[Hz] 



A 



S/N 
A 



is also defined Acp = cos 



cos Afl — COS 9y COS 9 
sin Qv sin 9 



Finally, 



the distance to the pulsar is D±. 

This formula was obtained in the form presented 
above first by Sago et al.jIS] in discussing the GW 
emission from gamma-ray bursts (GRBs). It is essentially 
equivalent to the formula obtained by Segalis and Ori 
[19( in an earlier paper in which they discussed the 
GW emission by ultrarelativistic sources. The main 
different in between is related to the detailed description 
of the angular distribution of the gravitational radiation 
emitted in each of the cases that those authors analyzed. 
The formula indicates that the GW strain depends on the 
Lorentz factor, pulsar speed and distance to the observer, 
and also on the strong beaming effect in the case of 
GRBs, as indicated by the factors (1 — cos 6), and AO, 
appearing in the denominator of the equation. These 
effects are no so much important in the case of Galactic 
RAPs, as is evident from the pulsar speed and angular 
size. 

This result, shown in FigJ^l states that the GW space- 
time perturbation is not strongly beamed in the forward 
direction ft, as opposed to the case of the electromagnetic 
radiation in gamma-ray bursts (also shown in the 
same Figure). Instead, the metric perturbation at the 
ultra-relativistic limit (not applicable to RAPs) has a 
directional dependence scaling as 1 + cos 9. In such a 
case, because of the strong beaming effect, as in gamma- 
ray bursts; for instance, the electromagnetic radiation 
emitted by the source over the same time interval is 
visible only inside the very small solid angle (6 ~ 7 -1 ) 2 , 
whereas the GW signal is observable within a wider 
solid angle; almost 2ir radians (see FiglJ). Besides, the 
observed GW frequency is Doppler blueshifted in the 
forward direction, and therefore the energy flux carried 
by the GWs is beamed in the forward direction, too. 
Quite noticeable, in the case of RAPs the GW signal will 
have its maximum for viewing angles 9 ~ tt/2, that is; for 
pulsar motions purely in the plane of the sky, as shown in 
Fig. This key physical property of the emission process 



J0014+4746 
J0024-7204C 
J0024-7204D 
J0024-7204E 
J0024-7204F 
J0024-7204G 
J0024-7204H 
J0024- 72041 
J0024-7204J 
J0024-7204N 
J0024-7204O 
J0024-7204U 
J0030+0451 
J0034-0534 
J0040+5716 
J0117+5914 
J0134-2937 
J0139+5814 
J0147+5922 
J0151-0635 
J0152-1637 
J0206-4028 
J0218+4232 
J0255-5304 
J0304+1932 
J0323+3944 
J0332+5434 
J0335+4555 
J0357+5236 
J0358+5413 
J0406+6138 
J0415+6954 
J0437-4715 
J0452-1759 
J0454+5543 
J0502+4654 
J0525+1115 
J0528+2200 
J0534+2200 
J0538+2817 
J0543+2329 
J0601-0527 
J0610-2100 
J0612+3721 
J0613-0200 
J0614+2229 
J0621+1002 
J0629+2415 
J0630-2834 
J0633+1746 
J0653+8051 



IT 



where 9 is the angle between v and n, that is the angleTABLE II: Characteristic amplitude and frequency, h c , / c (Hz), and 
between the unit spatial direction vector associated ttS/N = h c (f c )/h rma (f c ). 

velocity of the pulsar, V, and the unit spatial direction 

vector linking the source point to the observer location, Jname 

i. e. v-ri = cos 9). 9y is the particular (e. g. Vela _A 

Pulsar) viewing angle that a RAP velocity makes with 
the direction n to the observer. A9 is the angular 
diameter of the pulsar as seen by the observer. (For 
galactic distances it can be assumed nearly constant for 
each of the RAPs in the sample under analysis). Thus 
Aft ~ ir(A9) 2 is the solid angle over which the RAP is 
viewed. Besides, (3 — \V\/c, being \V\ the pulsar 3-D 
velocity V* n (t) defined above, and 7 the Lorentz factor; 
which can be taken here as ~ 1 because V+ n (t) <C c. It 



37.1239 7.99071x10" 
39.0641 1.81774xl0~ 24 
43.0236 1.13491 xl0~ 24 
38.1603 2.03054xl0~ 24 



40.5101 
40.1852 
39.0137 
38.98 
38.8167 



1.50882x10"^ 
1.57348xl0~ 24 
1.82934xl0 -24 
1.83713xl0~ 24 
1.87205 xl0~ 24 
37.9712 2.07818x10 
40.5101 1.50882x10 
38.1007 2.04539xl0~ 24 
79.288 2.32394xl0~ 25 
38.6853 9.31883xl0~ 24 
37.6624 2.31387xl0~ 24 
36.6742 5.64571 xl0~ 24 
36.5873 6.90089xl0~ 24 
38.4259 5.93842x10 
43.9134 2.67217x10 
38.4397 9.25843x10 
43.0987 6.85662x10 
38.5556 2.08228x10 
39.2696 1.45314x10 
39.0962 1.86966x10 
37.2766 1.14724x10 
37.1671 9.82703x10 
43.6786 4.89642x10 
53.156 1.38012x10 
37.143 3.12732x10 
52.7994 2.66806x10 
41.529 1.51989x10 
44.4139 3.13942xl0 -24 
42.4872 3.55097xl0~ 23 
38.6709 5.93204 xl0~ 24 
36.4967 1.65967xl0~ 23 
44.1655 2.817xl0~ 24 
40.4211 7.49499xl0~ 24 
36.9714 6.20109x10 
39.1178 4.33307x10 
39.3889 1.50976x10 
39.151 6.00026x10 
42.5181 5.41622x10 
38.8927 3.54204x10 
58.5628 1.34323x10 
71.1509 1.16838x10 
41.9176 1.30377x10 
67.466 4.79347x10 
37.6964 2.21027x10 
40.8252 1.17533x10 
40.4001 4.59112x10 
38.4132 5.25441x10 



-24 
24 



0.465967 
0.114416 
0.0810649 
0.123401 
0.100293 
0.103335 
0.114923 
0.115263 
0.116717 
0.125358 
0.100293 
0.124012 
0.0165996 
0.578055 
0.137876 
0.323257 
0.393722 
0.364667 
0.19087 
0.568847 
0.489759 
1.28516 
0.0921892 
1.1783 
0.673126 
0.57405 
0.349744 
0.09858 
0.182506 
0.190576 
1.04865 
0.224244 
2.53527 
0.367763 
0.943389 
0.201214 
0.496559 
0.359381 
0.273303 
0.962179 
0.378942 
0.386873 
0.221485 
0.0959447 
0.0834554 
0.091219 
0.0342391 
0.131882 
0.790392 
3.03935 
0.322503 
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TABLE II: results of full calculation of h c , fc and S/N. 



TABLE II: results of full calculation of h c , f c and S/N. 



Jname 


fc 


h c 


S/N 


Jname 


fc 


h c 


S/N 


A 


[Hz] 


A 


A 


A 


[Hz] 


A 


A 



J0659+1414 

J0700+6418 

J071 1-6830 

J0737-3039A 

J0737-3039B 

J0738-4042 

J0742-2822 

J0751+1807 

J0754+3231 

J0758-1528 

J0814+7429 

J0820-1350 

J0823+0159 

J0826+2637 

J0835-4510 

J0837+0610 

J0837-4135 

J0846-3533 

J0908-1739 

J0922+0638 

J0943+1631 

J0944-1354 

J0946+0951 

J0953+0755 

J1012+5307 

J1024-0719 

J1041-1942 

J1045-4509 

J1115+5030 

J1116-4122 

J1136+1551 

J1239+2453 

J1300+1240 

J1321+8323 

J1328-4357 

J1430-6623 

J1453-6413 

J1455-3330 

J1456-6843 

J1509+5531 

J1518+4904 

J2127+1155 

J1543-0620 

J1543+0929 

J1559-4438 

J1600-3053 

J1603-7202 

J1604-4909 

J1607-0032 

J1623-2631 

J1640+2224 

J1643-1224 

J1645-0317 

J1709-1640 

J1709+2313 

J1713+0747 

J1720-0212 



54.243 

71.8931 

42.4243 

80.003 

80.003 

41.1211 

37.6139 

74.4793 

38.44 

49.8372 

42.932 

41.3826 

65.7034 

36.7006 



37.0345 
39.0797 
42.1711 
38.8431 
40.8917 
36.5052 
50.0532 
36.5382 
64.7903 
53.4332 
38.7005 
36.6277 
37.5731 
39.0036 
45.2874 
42.1747 
40.3084 
35.5643 
36.5293 
41.7366 
38.6252 
37.0066 
46.1529 
45.4369 
40.5816 
67.6675 
42.3434 
44.2132 
43.8759 
42.6042 
46.2049 
57.8922 
40.8763 
73.4576 
38.4278 
52.6765 
36.5278 
39.5595 
74.3462 
45.5546 
66.3388 
41.9078 



9.24396x10"^ 
1.04344xl0~ 24 
5.59481 xl0~ 24 
1.37444xl0 -25 
1.37444xl0~ 25 
4.55418xl0~ 24 
8.28784xl0~ 24 
5.07004xl0~ 25 
2.40285xl0~ 24 
9.44338 xlO" 25 
1.27573xl0~ 23 
1.12452xl0~ 23 
7.55761 xl0~ 25 
3.31806xl0~ 23 
1.34572xl0~ 23 
1.56722xl0~ 23 
4.87755xl0~ 24 
2.1265xl0~ 23 
1.42168x10 
2.12817x10 
7.482x10 
5.0256x10" 
1.31884x10 
4.58178x10 
5.4258x10" 
1.71988x10 
4.21129x10 
3.23577x10 
1.62815x10 
1.73472x10 
8.4083x10 
2.84975x10 
2.88182x10 
1.66243xl0~ 23 
1.25324xl0~ 23 
1.02823xl0~ 23 
7.7727xl0~ 24 
5.94254xl0~ 24 
1.0231 xl0~ 23 
2.3723xl0~ 23 
1.25824xl0~ 24 
6.41853xl0~ 24 
4.34578xl0~ 24 
2.08026xl0~ 24 
3.56121 xl0~ 24 
1.64159xl0~ 24 
1.28057xl0~ 24 
7.19655xl0~ 24 
6.51078xl0" 25 
7.80249xl0~ 24 
2.48508xl0~ 24 
2.60936xl0~ 24 
7.80952xl0~ 24 
2.54331 xl0~ 25 
2.4967xl0~ 24 
9.09886xl0~ 25 
5.9315xl0~ 24 



23 
23 
-24 
24 
-23 
-24 
24 
-23 
-24 
-24 
-23 
-24 
-23 
23 
23 



0.660283 

0.0745312 

0.398562 

0.00981745 

0.00981745 

0.309597 

0.492893 

0.0362146 

0.147635 

0.0674527 

0.911233 

0.77176 

0.0539829 

1.90188 

0.961228 

0.910601 

0.307197 

1.50133 

0.887279 

1.43466 

0.425441 

0.358971 

0.750931 

0.32727 

0.387557 

1.06749 

0.240668 

0.192124 

1.02244 

0.123909 

5.9371 

1.88014 

1.57571 

0.946228 

0.871163 

0.636333 

0.451106 

0.424467 

0.730788 

1.58107 

0.0898742 

0.455935 

0.310413 

0.14859 

0.254372 

0.117257 

0.0914696 

0.484865 

0.0465056 

0.479169 

0.177506 

0.148511 

0.500944 

0.0181665 

0.178335 

0.0649919 

0.414854 



J1722-3207 

J1735-0724 

J 1738+0333 

J1740+1311 

J1741-3927 

J1744-1134 

J 1745-3040 

J1752-2806 

J1801-2451 

J1801-1417 

J1803-2137 

J1807-0847 

J1820-0427 

J1823+0550 

J1824-1945 

J1824-2452 

J1825-0935 

J1829-1751 

J1832-0827 

J1833-0827 

J1835-1106 

J 1836-0436 

J1836-1008 

J1840+5640 

J1841-0425 

J1844+1454 

J1850+1335 

J1857+0943 

J1900-2600 

J1902+0615 

J1905-0056 

J1906+0641 

J1907+4002 

J1909+1102 

J1909-3744 

J1910-5959A 

J1910-5959C 

J1911-1114 

J1913+1400 

J1913-0440 

J1915+1009 

J1915+1606 

J1916+0951 

J1917+1353 

J1919+0021 

J1921+2212 

J1926+1648 

J1932+1059 

J1935+1616 

J1937+2544 

J1939+2134 

J1941-2602 

J1944+0907 

J1946-2913 

J1946+1805 

J1948+3540 

J1952+3252 



41.6931 

41.8263 

52.7357 

42.0383 

40.304 

59.8904 

50.3561 

60.1363 

48.0157 

36.5256 

38.877 

45.9144 

37.5413 

48.5757 

40.9223 

42.3423 

49.139 

54.6043 

40.1473 

41.7718 

41.7953 

37.6347 

41.3733 

38.4105 

38. 

40.4852 
37.1645 
69.9194 
40.5195 
39.1075 
38.4463 
37.6858 
39.9779 
36.5478 
36.5355 
44.7141 
41.4165 
36.8491 
38.7526 
39.8321 
36.9264 
45.6686 
40.7829 
41.5323 
65.5118 
41.7992 
57.7089 
36.9246 
40.6174 
37.0187 
76.5964 
38.8781 
39.8356 
42.4317 
64.8701 
40.2965 
38.3669 



9.28191x10"^ 
6.49042 xl0~ 24 
1.49566xl0~ 24 
6.73415 xl0~ 24 
5.13915xl0~ 24 
3.77654xl0~ 24 
1.63708xl0~ 24 
1.16259xl0~ 24 
8.50308 xl0~ 25 
7.3598xl0~ 24 
5.109xl0~ 24 
1.24016xl0~ 24 
4.30552 xl0~ 24 
1.25897xl0~ 24 
2.43213xl0~ 23 
1. 19895 xl0~ 24 
3.6633xl0~ 24 
2.63599xl0~ 23 
5.06412x10 
7.73006x10 
1.43077x10 
3.43269x10" 
5.08732 xl0~ 24 
1.00787xl0~ 23 
3.26129xl0~ 24 
1.11881xl0~ 23 
4.66117xl0~ 24 
7.32693 xl0~ 25 
1.24926x10 
2.16284x10 
2.48948x10 
1.73977x10 
4.40831x10 
2.89161x10 
1.1301xl0~ 23 
1.20734xl0~ 24 
1.64278xl0~ 24 
7.30991 xl0~ 24 
3.7195xl0~ 24 
2.4511x10 
2.15708x10 
6.36111 xl0~ 25 
2.47226xl0~ 24 
1.59151xl0~ 24 
3.34207xl0~ 25 
9.49871 xl0~ 24 
3.11667xl0~ 24 
3.18859xl0~ 23 
3.16151xl0~ 24 
5.19088x10 
5.27885x10 
4.24401x10 
6.16328xl0~ 24 
8.50239xl0~ 24 
1.3907xl0~ 24 
3.1063x10" 
6.99324x10 



24 
24 
-23 
24 



23 
24 
-24 
24 
24 
24 



24 
-24 



24 
-26 
24 



24 
-24 



0.644202 

0.452622 

0.106833 

0.473193 

0.339001 

0.269753 

0.116934 

0.0830422 

0.0607363 

0.418842 

0.319274 

0.0885829 

0.255316 

0.0899267 

1.6414 

0.085163 

0.261664 

1.88285 

0.332107 

0.538018 

0.996669 

0.204318 

0.349026 

0.61854 

0.196949 

0.743005 

0.272255 

0.0523352 

0.830687 

0.136365 

0.152996 

0.103764 

0.28727 

0.16471 

0.643394 

0.0862383 

0.112883 

0.421542 

0.231327 

0.158855 

0.124784 

0.0454365 

0.165997 

0.109819 

0.0238719 

0.661766 

0.222619 

1.84443 

0.210985 

0.301413 

0.0037706 

0.26523 

0.399491 

0.60585 

0.099336 

0.204849 

0.428452 
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TABLE II: results of full calculation of he, fc and S/N. 



Jnamc 
A 



fc 

[Hz] 



A 



S/N 
A 



J1954+2923 

J1955+2908 

J1955+5059 

J1959+2048 

J2002+4050 

J2004+3137 

J2013+3845 

J2018+2839 

J2019+2425 

J2022+2854 

J2022+5154 

J2023+5037 

J2046-0421 

J2046+1540 

J2048-1616 

J2051-0827 

J2055+2209 

J2055+3630 

J2108+4441 

J2113+2754 

J2113+4644 

J2116+1414 

J2124-3358 

J2129+1210A 

J2129+1210B 

J2129+1210C 

J2129-5721 

J2145-0750 

J2149+6329 

J2150+5247 

J2157+4017 

J2219+4754 

J2225+6535 

J2229+6205 

J2229+2643 

J2235+1506 

J2257+5909 

J2305+3100 

J2308+5547 

J2313+4253 

J2317+1439 

J2321+6024 

J2322+2057 

J2326+6113 

J2330-2005 

J2337+6151 

J2354+6155 



45.8579 

43.7264 

40.8668 

36.7592 

41.0534 

40.3195 

35.0351 

67.6517 

43.7666 

38.4464 

41.1168 

36.6147 

36.5212 

37.6523 

38.8761 

67.0076 

57.4892 

41.9797 

44.3737 

39.1337 

38.4549 

37.9948 

53.7137 

36.6794 

42.2722 

37.31 

43.8083 

66.379 

40.4712 

36.9824 

40.3436 

39.1588 

41.5695 

37.6982 

41.5673 

43.6668 

41.3007 

39.1606 

37.1242 

43.4977 

51.393 

37.8303 

44.9277 

39.7623 

37.0468 

40.4429 

38.8978 



1.06684x10"^ 

9.59627xl0" 25 

1.40886xl0" 23 

9.00599xl0" 24 

3.17817xl0" 24 

2.74782xl0~ 24 

1.2161xl0~ 23 

9.09648xl0" 25 

5.66758xl0" 24 

6.63234xl0~ 24 

3.41183xl0~ 24 

6.76838xl0" 24 

3.44972xl0" 24 

4.05443 xlO" 24 

3.17922xl0" 23 

7.40902 xlO" 25 

1.00494xl0~ 24 

1.10321xl0~ 24 

9.36083xl0" 25 

1.56216xl0" 23 

3.59792xl0" 24 

3.80432xl0~ 24 

1.10886xl0" 23 

1.3541 xlO -24 

6.54473 xlO" 25 

1.08497xl0" 24 

2.01653xl0" 24 

2.02952xl0" 24 

4.24225xl0" 24 

2.51365xl0" 24 

4.42204 xlO" 24 

8.63068xl0" 24 

4.25672xl0" 23 

2.82253xl0" 24 

4.51013xl0~ 24 

4.51705xl0" 24 

4.20562xl0~ 24 

5.3719xl0" 24 

4.60111xl0~ 24 

5.53556xl0" 24 

1.6924xl0" 24 

5.16906xl0" 24 

6.10354xl0" 24 

4.82672xl0" 24 

2.29858xl0" 23 

2.95763xl0" 24 

6.1822xl0" 24 



0.762031 

0.0685448 

0.948881 

0.517451 

0.215521 

0.181364 

0.650212 

0.0649749 

0.404827 

0.407605 

0.231902 

0.386596 

0.196239 

0.241493 

1.9867 

0.0529216 

0.0717813 

0.0773579 

0.0668631 

0.985915 

0.221191 

0.229695 

0.792045 

0.0775484 

0.0463727 

0.0637445 

0.144038 

0.144966 

0.281582 

0.145743 

0.292128 

0.545228 

2.94121 

0.168426 

0.311606 

0.322647 

0.287776 

0.339384 

0.26831 

0.395397 

0.120886 

0.31007 

0.435967 

0.311996 

1.33621 

0.196108 

0.386652 



V. COMPUTING DETECTABILITY OF GW 
FROM RAPS 

One can discuss the detectability of the GW (burst) 
gnal of the rise fling of galactic RAPs by using the 
atching filter technique. This procedure is based on 



the definition of the signal-to-noise ratio (SNR), which is 
a quantity used as a criterium of detectability of a GW 
signal. The SNR depends on the features in the GW 
waveform of the process, the orientation of the source 
with respect to the GW interferometer (observer) and 
also on the source direction. It also depends on the total 
energy per unit frequency, dE/df(f) (the GW energy 
spectrum at the source), that is carried away from the 
source by the GW, and on the distance to the source D+. 
Here / is the average GW frequency for a given signal. 
Because of this angular dependence, the SNR is usually 
defined as 



SNR = (p 2 ) = 



dE 



PS h (/) df 



f df . (8) 



Symbol () refers to an average of SNR square for all 
the directions (angles) to the GW source. In other words, 
it is obtained from the average rms of the amplitudes of 
the GW signals for different orientation of the source and 
the interferometer. In our computations above in Section 
III we applied this definition to estimate the S/N ratio 
given in Table II. As the attentive reader may notice, in 
what follows we shall use an equivalent definition of the 
SNR to quantify the detectability of RAPs GW signals 
in the context of the formalism of Ref . [19| 

To complete this brief discussion on the SNR, we 
collect below the power-laws that according to most of 
the experimental studies with GW interferometers are 
what better describe the interferometer noise in each 
band of the spectrum. To this purpose , it is de fined the 
dimensionless quantity: h rrns (f) = yj fSh{f), wherein 
Sh(f) is the interferometer noise spectral density. Thus, 
the strain sensitivity for each frequency band is given by 



oo, 



(f) = \ h m (otf/fm) 3/2 , fs< f < fm/0i, 

S km, fm/a ^ f < afm, 

h m [II (a/m)] 3/2 , afm < f ■ 

(9) 

Understanding these laws is straightforward if one 
recalls that the interferometer noise curves depend 
basically on four parameters: 

• The threshold frequency f s : below of it the noise 
grows almost asymptotically. For detectors on 
Earth the noise source is the seismic gradients 

• The frequency f m , centered in the almost flat part 
of the noise spectrum 

• A critical amplitude h m , which is the minimum 
value of h rms (f) 

• And the dimensionless parameter a, which 
determines the width of the flat band of the noise 
curve 
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For instance, for Advanced LIGO the following set of 
parameters is estimated (see E. E. Flanagan and S. A. 
Hughes, [Phys. Rev. D 57, 8 (1998)]) 



LIGO Advanced 



= 10 Hz 
= 68 Hz 
1.6 

= 1.4 x 10" 



(10) 



->■■> 



Therefore, to know whether the RAPs GW signals 
can be effectively detected by the currently operative 
interferometric detectors LIGO, VIRGO, GEO-600, etc., 
one needs to compute the signal-to-noise ratio with 
respect to the strain sensitivity of these instruments. (We 
refer to this procedure as calibration). In order to do 
that we need first to compute the frequency at which the 
largest part of the GW emission is expected to take place. 
This is the so-called characteristic frequency / c , which is 
dependent on the detector design properties. 

Should the detector have a power spectral density 
Sh(f), have the characteristic frequency to be obtained 
from the frequency first moment of the time-dependent 
GW amplitude h(t), i. e., 



fc = 



s h (f) 



fdf 



Sh(f) 



-1 -1 



df 



(11) 



where (h(f)) is the Fourier transform of h(t) [32 |. 
and the symbol () stands for an average over randomly 
distributed angles of h(f), and it can be approximated 
as (\h(f)\) 2 = \h{f)\ 2 . 

The characteristic frequency f c is strictly related to 
the characteristic amplitude h c , which is the physical 
quantity associated to the GW pulse h(t) to be detected 
by a given instrument. It can be computed as 



h c (f c ) 



1 /2 

3 'o Wj ilHf)l2) 



(12) 



In this way we take h in the following form (see 
Ref.piGjl) 



\h\ 2 = ^^(l-cos(2nt m f)) 



(13) 



Here t m defines the rise time of the memory of the 
signal. This equation is obtained from the Eq.® given 
above. Thus, Eq. (ri~3f is the Fourier transform of the 
GW waveform, where Ah is given by (after average over 
angles in Eq.®) 



(14) 



with G the universal constant of gravitation, 7 Lorentz 
factor, M+ star mass, D+ distance to Earth, and (3 2 = 
\V\ 2 /c 2 . 

Therefore, the signal-to-noise ratio S/N is now 
computed using the relation 



S_ 

N 



hc(fc) 

Irms (.fc) 



(15) 



where h rms (f c ) = [f c S h (/c)] 1/2 is the average 
amplitude of the interferometer spectral noise at the 
characteristic frequency. 

Guided by Table ^TVl this allows us to compute the 
characteristic amplitude of the GW burst released during 
the short rise time of acceleration of each RAP for which 
those parameters are estimated. Since we are using 
only the transversal velocity of the pulsars, then we set 
6 ~ 90° in our calculations below. Our main results 
are presented in Figs. ( 4(a)] |4(b)| |5(a) ). Those Figures 
compare the GW signals from individual RAPs with the 
expected sensitivities of LIGO I, its projected advanced 
configuration LIGO II, VIRGO, GEO-600, and TAMA- 
300 observatories. 

At this stage, a note of caution to the attentive 
reader is needed. It regards with the appearance of 
the noise curves of TAMA, GEO and VIRGO GW 
observatories in figures from 4 through 9, in addition 
to the noise curves for LIGO I and LIGO Advanced, 
despite that the calibration of the expected signals 
in those figures (symbol: triangles, circles, etc.) is 
performed exclusively with respect to LIGO I and 
LIGO Advanced, independently. First, the reader must 
keep in mind that the position of those points in the 
diagram h c vs. f c will change when the calibration 
of the expect GW strain is performed with respect 
to the noise curve of a different detector, say TAMA, 
VIRGO, etc.. For instance, compare Fig. 4(a) and 
Fig. 5 (a) which are calibrated for LIGO Advanced and 
LIGO I, respectively. Second, thinking on the benefit 
of the Scientific Community operating observatories 
like TAMA, GEO and VIRGO, those noise curves are 
included having in mind the expectation that such 
detectors might gain more sensitivity in both frequency 
band and amplitude so as to be able to catch the GW 
signals produced by Milky Way RAPs. 

Finally, but not to the end, a direct interpretation of 
all the h c vs. f c diagrams in figures from 4 through 9 
is that, exception done for the TAMA-300 GW detector, 
the signals (points) plotted in those figures appear to 
be not very much sensitive to the noise curves of the 
instruments LIGO I, LIGO Advanced, VIRGO and GEO- 
600. This feature justifies to include those noise curves 
in such figures. 
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(a) LIGO I, LIGO Advanced, TAMA, VIRGO and (b) Idem as Fig.( |4(a)| , but now for parameters a = 10 6 km 

GEO600 strain sensitivities and the h c and f c s -2 and V^ick 10* km B —1 , calibrated with respect to 

characteristics of the GW signal produced by each of the Advanced LIGO. It is possible to see that these signals have 

212 pulsars calibrated with respect to LIGO Advanced, for high enough chance to be detected, 
parameters a = 10 4 km s~ 2 and Vkick ~ 10 4 km s — 1 . 



FIG. 4: Color Online. Pulsar GW signals 




different parameters a and Vkick- ( a )> (b) 

VI. MORE REALISTIC ANALYSIS OF GW 
EMISSION FROM RAPS 

We can start with by recalling both the inertia 
(which is large for a neutron star progenitor) and linear 
momentum conservation laws, which should dominate 
the physics of the SN explosion that gives the kick to 
the neutron star. Thus, one can write 



M S nV sn = M NS V NS ~ MrsnVrsn , (16) 



10 



•24| 

10 



100 
f [Hz] 



1000 



(a) Idem as Fig. | |4(a)| l, but now for parameters a = 10 6 km 
s — 2 and Vkick ~ 10* km s — 1 , and calibrated against LIGO 
I. Clearly also these signals might be detected 



FIG. 5: Color Online, 
parameters a and V. (c) 



Pulsar GW signals for different 



where Msn, Mns, and Mrsn are the SN, neutron 
star, and ejected material mass, respectively. Be aware 
of the role of the "— " signal in front of the term for 
the ejected material. It is responsible for a correct 
description of the physics at the explosion. Without it, 
one would obtain immaginary values for the searched 
parameter Vns- Now, by assuming parameters of a 
typical progenitor star Msn — 10 M Q , Vsn = 15 km 
s _1 , Mjvs = 1.4 M Q for a canonical neutron star, and 
for the material in the ejected envelope Mrsn = 8.6M Q , 
and expansion velocity Vrsn = 3000 km s _1 (see Table- 
IIV[) . one obtains that the kick at birth velocity reads: 
Vo = Vkick ~ 10 4 km s -1 (compare to velocities taken 
from ATNF in Fig 1(c) I. The discrepancy between them 
will be briefly discussed below. 




(a) Efficiency case e = 0.40 for parameters a = 10 6 km (b) Efficiency case t = 0.10 for parameters a = 10 6 km 

s — 2 and Vkick ~ 10 4 km s_ 1 i calibrated against LIGO s~ 2 and Vkick ~ 10 4 km s — 1 , calibrated against LIGO 

Advanced. Advanced. 
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1000 



f [Hz] 

c 



(c) Efficiency case e = 0.05 for parameters a = 10 6 km 
s~ 2 and Vkick ~ 10 4 km s —1 , calibrated against LIGO 
Advanced. 
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(d) Efficiency case e = 0.01 for parameters a = 10 6 km 
s~ 2 and Vkick ~ 10 4 km s _1 , calibrated against LIGO 
Advanced. 



FIG. 6: Color Online. Pulsar GW signals for parameters: a = 10 4 km s 2 and Vkick ~ 10 4 km s 1 



A. Case a: Centered kick V = Vkick and e — 1 

We first consider that the whole energy obtained from 
the conservation law discussed above goes into the pulsar 
kick, imparting to it no rotation at all. This is a very 
idealistic situation. Nonetheless, it gives us an idea on 
the potential detectability of those GW signals from the 
flinging of a pulsar out of the SN cocoon. In other words, 



in this first approach we assume that the efficiency, e, 
that relates the effective RAP velocity V (which enters 
the calculation of h c ) and the kick speed Vkick is e = 1 
(In the other cases analyzed below these velocities will 
be related through the expression: V = e Vkick- 
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TABLE III: ATNF sample of pulsars having 2-Dim speed 
greater than 10 3 km s _1 . The current high velocities suggest 
that a much higher speed was imparted at the kick at birth. 



1 f i > \ iy i i_ \ 








[kpcj 


[K1T1 S 


O yjO-^O] _L _L _L tj 


7.68 


1 1 01 Q8*l 
llul. iJOO 


J1509+5531 


2.41 


1104.013 


J1824-1945 


5.20 


2482.989 


J1829-1751 


5.49 


3945.954 


J2013+3845 


13.07 


2521.130 


J2149+6329 


13.65 


1113.385 


J2225+6535 


2.00 


1729.770 


B2011+38 


???? 


1624.0 


B2224+65 


???? 


1608.0 


B1830-1945 


5.20 


740.0 


*RX J0822-4300 


2.0 


1600.0 



*Ref. P. F. Winkler & R. Petre, Direct measurement of neutron 
star recoil in the Oxygen-rich supernova remnant Puppis A, Ap.J. 
670, 635 (2007) 



supposed to be expended in making the pulsar to spin. 
This parametrization can be translated into a power- 
law relationship between the GW amplitude, h, and the 
efficiency, e, of energy conversion into pulsar kick, which 
has the form: 



hoce 2 . (17) 

This relation stems from the quadratic dependence of 
the GW amplitude with the (3 parameter appearing in 
Eq.([7]) above. In this way, one can recover the whole 
spectrum of RAPs velocities as observed today and listed 
in ATNF catalog. To exemplify, we compute signals 
for values of e = 0.01,0.05,0.10, and 0.40, which can 
reproduce the set of velocities I4ick = 100 km s _1 , Vkick = 
500 km s" 1 , Vk ick = 10 3 km s _ \ and Vudk = 4 x 10 3 km 



s , as shown in Fig.-( 1(c) ). For this case, the results are 
presented in Figs.© UTM M ■ To the end of the 

paper 

a brief discussion on the pulsar periods provided by the 
ATNF catalog is given. 



1, Fling with a = 10 4 km s 2 , Vkick = 10 4 km s 1 

Based on the dynamical analysis presented in the 
Section- Hill we suppose here that the short rise time kick 
imparts to the pulsar an acceleration a = 10 4 km s~ 2 , 
and a birth velocity Vkick = 10 4 km s _1 . The GW signal 



obtained for these parameters is presented in Figs. (4(a) 



2. Fling with a = 10 6 km s -2 , Vki ck = 10 4 km s" 1 

Now one can take a step further and consider that the 
birth acceleration is actually more larger than the one 
currently inferred from the RAPs galactic parameters. 
Several physical mechanisms can act together to impart 
an acceleration as large as a — 10 6 km s~ 2 , see for 
instance Refs.(23l, l24l [26| . and references therein. Such 
a large thrust can be estimated as a ~ I4i C fe/AT mec ft, 
where (10 _1 < AT mec h < 10~ 3 ) s, depending on the 
specific mechanism driving the fling. The GW signal 
obtained for these parameters is presented in Fig. (5(a) ). 



1. Efficiency e = 0.05 

The following figures illustrate the results obtained 
for this efficiency: Fig. (7(c) I with a — 10 6 km s~ 2 and 



Fig. (6(c) ) with a = 10 4 km s 2 for advanced. Fig. (8(c)) 
with a = 10 6 km s~ 2 and Fig.( |9(c)[ ) with a = 10 4 kmP 
para LIGO I. 



2. Efficiency e = 0.10 

The following figures illustrate the results obtained for 
this efficiency: Fig. (7(b)) with a = 10 6 km s~ 2 and 
Fig. ([6(b)]) with a = lCFkm s~ 2 for Adva nced L IGO. 
Fig. |9(EJ ) with a = 10 6 km s~ 2 and Fig. ( [8(b)) with 
a = lCFkm s~ 2 for LIGO I. 



3. Efficiency e = 0.40 

The following figures illustrate the results obtained for 
this efficiency: Fig. (7(a)) with a = 10 6 km s~ 2 and 
Fig.(|6(a)| with a = RT km s~ 2 for Adva nced L IGO. 
Fig.(|9(a)f with a = 10 6 km s~ 2 and Fig. ([8(a)]) with 



km s 
a = 10 6 km s~ 2 
a = lCFkm s~ 2 for LIGO I. 



B. Case b: Off-centered kick 



Finally, we analyze the most astrophysically realistic 
kicks at birth: a situation where the pulsar receives 
both a translational and a rotationl kick. In this 
case, the energy transferred to the translational kick 
can be parameterized by an efficiency e, which ranges 
from to 1. Of course, this efficiency represents our 
difficulty in explaining how a specific kick mechanism 
becomes the kick driver. The remaining energy is 



VII. DISCUSSION AND CONCLUSIONS 

A. Computed rotation periods vs. Periods of 
ATNF sample of RAPs with V = 400 - 500 km s _1 

To the end, in order to make our analysis self- 
consistent we next discuss the relation between the 
translational speed and rotation frequency of ATNF 
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TABLE IV: High expansion speeds of Call and Sill lines of the ejected material in some supernovae. Data taken from A. 
Balastegui, P. Ruiz-Lapuente, J. Mendez, G. Altavilla, M. Irwin, K. Schahmaneche, C. Balland, R. Pain and N. Walton. arXiv: 



astro-ph/0502398^1 21 Feb 2005 



Event Call Expansion Speed Sill Expansion Speed Days after maximum 
[km s _1 ] [km s _1 ] [days] 



SN 20021i 
SN 20021j 
SN 20021p 
SN 20021q 
SN 20021r 
SN 20021k 



17200± 500 
12300± 300 
13800± 900 
18900± 1400 
12800± 300 



8900± 300 
10400± 400 

8600± 400 
14500± 500 



-7 
7 
3 
-7 
10 
-5 



pulsars. To do this, we select, as an example, the 3- 
D average pulsar velocity as determined by Hobbs et 
al. (see Table-TV]). On this basis we determine 

the theoretical rotation periods that one could expect 
to find for that velocity sample if is factual that most 
of the explosion energy, as inferred from the momentum 
conservation law in Eq. (|16|) . goes into the pulsar rotation 
as discussed above. (Of course we compare with the 
actual periods measured by the ATNF surveys). Hence 
the torque applied to the just-born neutron star reads 



V lM NS R 



2E k 



H 



(18) 



where E k is the kinetic energy of each velocity chosen 
for the present analysis, V pulsar velocity, Mats = 1.4 
M©, Rats = 20 km and d; = 25 km are the neutron 
star mass, radius and torque lever. Thus, for Kns = 20 
km and velocity V — 350 km s _1 we obtain p= 0.0182 
s, while for V = 433 km s -1 one gets p= 0.0182 s. 
Meanwhile, for Kns — 30 km and V = 350 km s _1 we 
have p= 0.0411 s, while for V = 433 km s _1 one gets 
p= 0.0332 s. 

One can then compare these results with those 
presented in Table-fVl It appears that the theoretical 
prediction is about an order of magnitude smaller (much 
shorter periods) than the observed, and measured by the 
ATNF surveys. This discrepancy may suggest that 

• 1) perhaps a large part of the initial rotational 
energy is dissipated by processes as differential 
rotation or magnetic braking, especially if processes 
as the magneto-rotational instability [27| indeed 
drives a large class of supernova explosions [22j ■ 

• 2) it may also happen that convection and neutron 
fingers also becomes highly effective dissipative 
processes, thus reducing the pulsar initial spin to 



the level that we are currently finding in pulsar 
surveys. 

In summary, from astronomical statistics the number 
of observed Galactic RAPs is definitely much larger than 
the corresponding to relativistic NS-NS binaries. The 
typical lifetime of canonical pulsars is ~ 10 Myr, which 
is defined as the duration of the radio emitting phase 
for objects with B ~ 10 12 G, or the time needed for 
the pulsar to emit less radiation or being entirely turned 
off. Hence, with a statistical Galactic rate of SNe, and 
the likely fling out rate of just-born pulsars, of about 1 
per every 30-300 years, one expects that the GW signal 
of the kick at birth be observable as aftermath of the 
SN core collapse. In other words, the RAP GW signal 
may become a benchmark for looking for GW signals 
from the SN gravitational implosion. Besides, if the GW 
signal from the kick to a RAP were detected, for instance, 
from a source inside a globular, or from either an open 
star cluster [T(| , or a binary system harboring a compact 
star, such an observation might turn these RAPs the 
most compelling evidence for the existence of Einstein's 
gravitational waves. 

We conclude by stating that the detection with 
advanced LIGO-type interferometers of individual GW 
signals from the short rise fling of Galactic RAPs is 
indeed possible. 
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